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ABSTRACT 


The  Free  Electron  Laser  (FEL)  is  a  potential  solution  for  the  U.S.  Navy's  anti-ship 
missile  point  defense  by  providing  an  evolutionary  increase  in  Aveapon  accuracy.  To  become 
an  effective  weapon,  the  FEL  will  need  to  provide  an  average  optical  power  of  approximately 
one  MW.  Towards  this  goal,  the  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  in 
Newport  News,  Virginia  is  constructing  the  first  kW  FEL,  and  desires  to  improve  the  design 
to  20  kW  while  maintaining  less  than  6%  energy  spread.  Using  a  klystron  undulator  is  one 
potential  way  to  accomplish  this.  Given  design  parameters  of  a  proposed  free  electron  laser 
by  TJNAF,  this  study  quantifies  via  simulation  the  behaviors  of  gain,  power  and  energy 
spread  as  functions  of  desynchronism  and  a  klystron's  dispersive  strength.  Specifically,  it 
shows  that  a  conventional  imdulator  appears  capable  of  meeting  all  TJNAF  design 
requirements. 
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I.  INTRODUCTION 

A.  THE  DEFICIENCY  IN  CURRENT  ANTI-SHIP  MISSILE  POINT  DEFENSE 

1.  Current  Point  Defense  System  Description 

Since  1978,  the  Phalanx  Close-In  Weapon  System  (CIWS)  has  been  the  U.S. 
Navy's  last  line  of  defense  against  anti-ship  cruise  missiles.  CIWS  is  currently  installed 
aboard  all  combatants,  amphibious  landing  and  logistics  support  ships.  CIWS  uses  a  20- 
mm  Gattling  gun  with  a  microwave  radar  that  tracks  both  the  incoming  missile  as  well  as 
the  outgoing  CIWS  rounds.  This  allows  the  moimt  to  autonomously  adjust  the  stream  of 
rounds  onto  the  target.  Engagements  can  begin  out  to  2  km  with  a  minimum  engagement 
range  of  100  m. 

2.  Performance  Appraisal  Using  Simulation 

Just  how  effective  is  this  system?  CIWS  has  never  been  challenged  in  combat. 
We  must  therefore  resort  to  simulation  to  predict  its  performance  against  real  world 
threats.  Modeling  the  path  of  CIWS  rounds  is  a  simple  problem  in  classical  physics.  For 
the  purpose  of  this  analysis,  the  effects  of  gravity  and  air  drag  are  included,  and  perfect 
tracking  by  CIWS'  radar  is  assumed. 

Every  gun  system  suffers  from  an  effect  called  dispersion.  It  is  the  error  incurred 
when  rounds  do  not  go  precisely  where  they  are  directed.  The  reasons  for  this  effect  are 
numerous:  barrel  vibrations,  variations  in  shell  manufacturing,  etc.  Despite  reduction 
efforts  in  this  area,  CIWS  still  experiences  a  dispersion  of  3  millirads  or  0.17  degrees  [1]. 
It  may  not  soimd  like  much,  but  it  has  a  large  affect  upon  performance. 
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A  typical  CIWS  burst  is  approximately  200  rounds.  Given  a  stationary  missile  with  a  20 
cm  radius  (typical  of  today's  cruise  missile  inventory)  at  750  m  range  (less  than  half  the 
maximum  engagement  range),  only  3  rounds  out  of  200  actually  hit  [1]!  Figure  1 
graphically  shows  these  results. 


CIWS  Dispersion 

CIWS  Height  25m,  Dispersion  3mrads,  Missile  Height  5m 
Missile  Range  750m,  200  Rounds,  Seed=28 


Figure  1.  CIWS  Dispersion 


This  makes  the  intercept  probability  of  any  one  round  at  750  m  just  1.5%.  Figure  2 
shows  the  intercept  probability  behavior  as  a  function  of  decreasing  range  [1].  The 
maximum  probability  is  just  over  35%  at  100  m  (CIWS'  minimum  engagement  range). 
Beyond  700  m,  the  probability  of  hit  drops  off  to  almost  nothing. 
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Single  Round  Intercept  Probability 
CIWS  Height  25m,  Dispersion  3mrads,  Missile  Height  5m, 
1 0000  Rounds  Fired  Every  1 00m,  Seed=28 


Figure  2.  Single  Round  Intercept  Probability 

The  previous  analysis  has  focused  on  hitting  the  missile  once.  Unfortunately, 
missiles  rarely  die  from  the  first  hit.  Rather,  it  takes  an  average  of  six  to  ten  hits  to  either 
destroy  the  missile’s  aerodynamics  or  ignite  its  warhead  [1].  To  continue  with  our 
analysis,  let  us  assume  eight  hits,  and  that  the  missile  is  moving  at  300  m/s  (slow  by 
today’s  standard)  on  a  constant  trajectory  inboimd  to  the  ship.  At  what  range  is  this 
missile  typically  destroyed? 

The  CIWS  firing  rate  is  50  rounds/s,  and  each  round  has  an  initial  velocity  of 
1,200  m/s.  This  means  that  an  incoming  missile  encounters  approximately  21  roimds 
with  every  100  m  section  in  decreasing  range  to  the  ship.  The  number  of  CIWS  hits  on 
the  missile  at  a  particular  range  is  equal  to  21  rounds  multiplied  by  the  siiigle  round 
intercept  probabilities  shown  in  Figure  2.  The  accumulated  hit  count  as  the  missile 
approaches  the  ship  reaches  eight  at  approximately  150  m.  If  we  simulate  debris 
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trajectories  from  a  typical  warhead  at  a  height  5  m  above  the  water  and  150  m  range, 
more  than  half  of  the  debris  strikes  the  ship  [1]! 

So  for  even  a  quite  simple  threat,  the  ship  will  suffer  damage.  In  today's  warfare 
of  zero  tolerable  losses,  any  amount  of  ship  damage  could  force  its  withdrawal  from  the 
theater  of  operations.  Given  today’s  sophisticated,  anti-ship  cruise  missile  environment, 
CIWS  can  not  be  counted  upon  to  capably  defend  our  fleet.  Therefore,  we  should  pursue 
research  on  a  point  air  defense  system  that  is  much  more  accurate  with  longer  range. 

B.  POTENTIAL  ANSWER:  THE  FREE  ELECTRON  LASER 

1.  Basic  Description 

A  potential  long-term  answer  to  the  U.S.  Navy’s  point  air  defense  problem  is  a 
high-energy  laser.  One  candidate  is  the  Free  Electron  Laser  (FEL).  In  its  simplest 
description,  a  FEL  is  made  of  two  major  components:  an  electron  accelerator  and  an 
undulator.  The  accelerator  propels  electrons  to  near  the  speed  of  light.  The  relativistic 
electrons  then  enter  the  undulator,  which  consists  of  a  series  of  alternating  magnets.  If 
we  assume  a  planar  undulator  (2D  only),  these  alternating  magnets  produce  a  sinusoidal 
magnetic  field  along  the  undulator  axis.  As  the  electrons  pass  by  each  magnet,  they  are 
deflected  slightly  from  side  to  side.  This  periodic  action  converts  the  straight  electron 
beam  path  into  an  oscillating  one.  As  the  electrons  oscillate  in  the  presence  of  light,  they 
bunch  together  and  emit  coherent  light  radiation. 

There  are  two  types  of  free  electron  lasers:  oscillators  and  amplifiers.  In  an 
oscillator  (Figure  3),  the  optical  field  (seen  in  gray)  is  stored  between  two  mirrors  with 
the  undulator  in  the  middle.  With  each  pass  of  the  optical  wave  through  the  undulator,  it 
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receives  more  and  more  energy  from  the  electron  beam.  Between  one  and  ten  percent  is 
extracted  from  the  optical  cavity  to  make  a  weapon! 


Electron  Beam 


Undulator 


Electron  Dump 


Figures.  Oscillating FEL 

An  amplifier  is  similar,  but  lacks  the  mirrors  at  either  end  of  the  imdulator.  Here, 
the  optical  wave  makes  only  a  single  pass  down  the  undulator.  This  design  works  well  if 
the  optical  field  already  exists,  the  FEL  gain  is  large,  and  the  output  power  is  enough  to 
destroy  any  mirrors. 


2.  Potential  Advantages  Over  The  Current  Point  Defense  System 

While  FEL  offers  several  advantages  over  CIWS,  the  primary  ones  are  increased 
range  and  improved  accuracy.  Although  CIWS  has  a  maximum  engagement  range  of  2 
km,  the  previous  analysis  shows  that  consummation  of  even  a  simple  engagement  may 
not  occur  imtil  a  few  hundred  meters.  In  addition  to  subjecting  the  ship  to  missile  debris 
at  this  short  kill  range,  CIWS  is  unlikely  to  have  sufficient  time  to  engage  a  second 
trailing  missile.  On  the  other  hand,  the  FEL's  effective  range  is  from  5  to  10  km.  This  is 
a  tremendous  tactical  advantage  since  it  prevents  ship  damage  from  debris  and  greatly 
improves  our  ability  to  engage  a  second  and  even  third  successive  threat. 


5 


A  second  advantage  is  improved  accuracy.  Laser  weapons  do  not  suffer  from  dispersion. 
Through  the  use  of  optics,  we  can  focus  a  laser's  energy  within  the  cross-sectional  area  of 
a  missile  even  at  5  to  10  km. 

3.  Potential  Advantages  Over  Chemical  Lasers 

The  U.S.  Army  and  Air  Force  are  developing  high-energy  chemical  lasers  for 
battlefield  air  defense  and  Theater  Ballistic  Missile  Defense  (TBMD)  respectfully.  Why 
not  choose  such  a  laser  for  maritime  defense?  There  are  two  main  disadvantages.  First, 
the  chemical  process  chosen  for  each  laser  application  dictates  its  wavelength.  This 
makes  the  weapon  inflexible  to  changes  during  and  after  the  design  process.  A  second 
disadvantage  is  that  high  power  chemical  lasers  produce  large  quantities  of  hazardous 
waste  that  can  not  be  disposed  of  at  sea.  This  greatly  complicates  tailoring  such  a  laser 
for  maritime  employment.  For  these  reasons,  chemical  lasers  may  not  be  the  Navy's 
answer. 

A  free  electron  laser  does  not  produce  hazardous  waste,  and  its  wavelength  is 
tunable  by  three  means:  changing  the  imdulator's  wavelength  (distance  between 
magnets),  the  imdulator's  magnetic  field  amplitude  or  the  accelerator's  electron  beam 
energy.  This  provides  great  flexibility  in  optimizing  a  laser  within  its  operating 
environment. 

4.  Power  Required  to  Destroy  a  Missile 

Just  how  much  optical  power  is  required  to  destroy  a  missile?  For  ship  self- 
defense,  the  goal  is  to  cause  structural  failure  of  the  missile's  nose  cone  rather  than  ignite 
its  warhead.  Let  us  assume  that  the  nose  of  the  missile  is  approximately  1  L  in  volume 
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consisting  of  about  3x10^^  atoms.  In  atypical  solid,  we  need  approximately  1  eV 
(electron-volt)  of  energy  to  remove  each  atom.  This  translates  to  5  MJ  of  total  energy.  If 
we  permit  2  to  3  seconds  per  engagement,  this  means  we  need  1.5  to  2  MW  of  average 
power  delivered  at  the  missile  [1]. 

5.  Present  Stage  of  Research  &  Development 

How  close  are  we  to  achieving  a  MW  PEL  suitable  for  ^pboard  deployment?  In 
the  1980's,  the  Strategic  Defense  Initiative  (SDI)  goal  was  to  design  a  10  MW  PEL  with 
no  intermediate  steps  at  smaller  powers.  Scientists  completed  the  design,  and  if  built  it 
would  have  been  quite  large  and  costly.  With  a  change  in  national  priorities,  the  U.S. 
Government  cancelled  SDI  before  completion  of  the  PEL  project.  Since  then,  several 
laboratories  have  worked  with  operational  PELs  in  the  10  watt  range.  Currently,  the 
Thomas  Jefferson  National  Accelerator  Pacility  (TJNAP)  in  Newport  News,  Virginia  has 
undertaken  the  next  step  by  constructing  a  kW  PEL.  Operation  begins  in  the  summer  of 
1998. 

C.  PURPOSE  OF  THIS  THESIS 

Once  TJNAP  has  successfully  proven  the  kW  PEL,  their  next  design  step  is  to 
increase  output  to  20kW.  There  are  several  design  initiatives  worth  exploring  to  achieve 
this.  One  is  to  use  a  spedal  type  of  undulator  called  a  klystron.  The  purpose  of  this 
thesis  is  to  quantify  the  steady  state  gain  and  final  power  behaviors  of  a  proposed  20  kW 
klystron. 
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11.  FEL  THEORY 


A.  ELECTRON  BEHAVIOR  WITHIN  THE  UNDULATOR 

To  understand  how  a  FEL  works,  it  is  first  important  to  comprehend  the  electron 
behavior  within  the  undulator  [2],  They  travel  at  relativistic  velocities,  and  their  motion 
is  described  using  the  principles  of  classical  physics.  For  the  undulator,  let  us  use  the 
coordinate  system  shown  in  Figure  4: 

X 

z 

Figure  4.  Undulator  Coordinate  System 

The  imdulator's  circularly  polarized  magnetic  field  can  be  described  by  the 
equation: 

Bq  =  5[cos  k^z,  sin  k^z,  0]  (1) 

where  B  is  the  magnetic  field  amplitude.  kQ=lnl  Xq  is  the  undulator  wave  number,. 

Xq  is  the  undulator  wavelength  (distance  between  magnets),  and  z  is  the  distance  down 

the  undulator.  The  generated  optical  field  within  the  undulator  is  a  circularly  polarized 
plane-wave  composed  of  its  own  electric  and  magnetic  fields.  These  optically  induced 
fields  can  be  described  by  the  equations: 
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£  =  £[cos'P-sin'F,0] 


and  B  =  £'[sin  T,  cos  T,0] 


(2) 


E  is  the  electric  and  magnetic  field  magnitudes  while  =  kz -cot +  ^ .  k  =  co  I  c'lsiht 

optical  field  wave  number,  co  is  the  optical  wave's  angular  frequency,  and  <f>  is  the  optical 
wave's  phase.  The  speed  of  light  in  a  vacuum  is  c,  and  t  is  the  time. 

The  Lorentz  force  equation  governs  electron  motion  in  the  presence  of  electric 
and  magnetic  fields, 

f  =  and  =  (3) 

where  p  =  ytm  is  the  electron  momentum,  m  is  the  electron  mass,  e  is  the  electron 
charge,  v  is  the  electron  velocity,  >0  =  v  /  c  is  the  dimensionless  electron  velocity  and 
y  =  {^-  the  Lorentz  factor.  From  Einstein’s  relativistic  theory,  the  electron 

energy  is  ymc^ .  Finally,  +  5  is  the  total  magnetic  field  within  the  undulator. 

Substituting  the  definition  for  momentum  into  the  Lorentz  force  equation,  we  can  get  (3) 
into  the  form: 

M  =  +  and  (4) 

at  me  at  me 

We  now  insert  the  electric  and  magnetic  fields  of  (1)  and  (2)  into  (4),  and  simplify 

the  results  by  splitting  P  into  two  components:  P^  Pz  Pz\^^ 

arrive  at  three  components: 

iMA  =  _  (£(i  _  )[cos  T,-  sin  4^,0]  +  Bp^  [-  sin  z,  cos  z,0]) ,  (5) 

dt  me 
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=  cos'F-y^,  +  smkoZ-p^  cosk^z],  (6) 

CIT  TfYtC 


and  ^  =  — ^My^v.  cosT-y^y  sin'F). 
at  me 


(7) 


To  understand  electron  behavior  within  the  undulator,  we  need  only  use  equations 
(5)  and  (7)  because  (6)  can  be  derived  from  (5)  and  (7).  Concerning  equation  (5),  we  can 
make  several  simplifying  approximations.  First,  » 1  so  £(1  -  y5y  )  «  when 
Ef=iB.  This  allows  us  to  ignore  the  first  term  in  (5)  so  that 

eBp 


■  [-  sin(^^  z,  cos(^o  z),0] . 


dt  me 

After  integrating  (8),  the  component  equations  become: 


Px  ,  cosikoz)  and  ygy  ,  sin(A:oz). 


(8) 


k^yme^ 


k^yme 


(9) 


To  further  simplify  (9),  the  dimensionless  imdulator  parameter  is  AT  =  eBIkome^  so: 


Px  =-{EI y)cos{kQz)  and  Py  =-{K ! y)sm(kQz)  (10) 

The  next  step  is  to  now  insert  (10)  into  (7)  using  the  trigonometric  identity: 
cos  A  cos  5  -  sin  ^  sin  B  =  cos(^  +  B) , 


^  =— cos(^oZ  +  'F), 
at  ymc 


(11) 


where  k^z  +  ^i*  =  kQZ  +  kz-a)t  +  ^  =  (kQ +k)z-ca  +  ^ .  Define  ^  =  (,ky  +ko)z-cot  as 
the  electron  phase  in  the  combined  imdulater  magnetic  field  and  the  associated  optical 
fields  so  that  (1 1)  becomes: 

dy  .  eKE 
—  =  r  = - cos{C  +  <P). 

dt  yme 
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(12) 


Equation  (12)  looks  similar  to  a  simple  pendulum  equation,  but  to  get  it  into  the  right 
form,  we  need  to  express  y  in  terms  of  ^ .  Recall  the  definition  for  the  Lorentz  factor: 
l/y"  =1-J3^  From  (9),  / y"  so  that: 


Since  ^  =  (k^,  +  kp)z  -  cot , 


^  =  {ko+k)cp2.  (14) 

When  (14)  is  substituted  into  (13),  the  result  is: 

f  _  (  7^ 

y  (Ato +A:)c|^1  +  A:^ 

To  simplify  any  further,  we  need  to  shift  our  focus  to  the  interaction  of  electrons 
and  photons.  Photons  overtake  the  electrons  as  they  both  proceed  down  the  undulator. 
When  one  wavelength  of  light  X  overtakes  an  electron  just  as  the  electron  traverses  one 
undulator  wavelength  ,  the  FEL  is  at  resonance.  Finding  how  the  imdulator  and 
optical  field  wavelengths  are  related  at  resonance  will  enable  us  to  simplify  (15),  which 
will  in  turn  allow  us  to  write  (12)  in  the  desired  simple  pendulum  form.  The  time 
required  for  an  electron  to  traverse  one  undulator  period  is  Iv^  =Xq  I  2^ .  In 

the  same  time,  the  traveled  distance  of  a  photon  is  A  =  ct^  ~  fiz)  ■ 

Substituting  the  first  condition  into  the  second  yields  the  relationship: 
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(17) 


2X  \  +  K^ 

Reverting  back  now  to  (15),  we  can  simplify  it  by  substituting  in  equation  (17), 


y 


Cr 

(  k  ^ 

{k^  +  k)c 

(Icq  +  k)c 

(18) 


The  undulator  wavelength  (cm)  is  much  greater  than  the  laser's  optical  wavelength  (pm) 
so  k»ko.  This  allows  (1 8)  to  simplify  to: 


IkpC 

Finally,  inserting  (19)  into  (12),  the  pendulum  equation  of  electron  motion  within  the 
undulator  is: 


(19) 


IknsKE 

C  =  - cos(C  +  ^).  (20) 

y  m 

It's  helpful  when  designing  a  FEL  to  use  dimensionless  parameters  when  possible. 
Such  is  the  case  when  dealing  with  time.  For  an  average  electron,  the  time  it  spends  in 
the  undulator  is  t  =  Lfv2  =LI  ^  where  L  is  the  undulator  length  and  » 1 . 

The  dimensionless  time  is  defined  as  x  -  ct  /  L  so  that  t=0  at  the  undulator  beginning  and 
T=1  at  the  end.  Rewriting  (20)  in  terms  of  the  second  derivative  with  respect  to 
dimensionless  time  yields; 


^  =  \a\cos{C+<P),  (21) 

where  the  dimensionless  optical  field  is  |a|  =  4nNeKLE I y^mc^  ,ihs.  number  of  imdulator 
periods  is  N  =  LI Xq,  and  the  open  dot  represents  the  derivative  with  respect  to  t . 
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B.  OPTICAL  WAVE  BEHAVIOR  WITHIN  THE  UNDULATOR 


The  previous  section  has  detailed  electron  behavior  within  the  undulator.  Now, 
attention  is  turned  to  how  the  optical  wave  develops  within  the  undulator  [3],  As  stated 
earlier,  the  optical  field  is  a  circularly  polarized  plane-wave.  The  wave's  vector  potential 
is: 

A  =  [sin  'F,  cos  ^,0].  (22) 

E(z,t)  is  the  optical  wave's  electric  field  magnitude  and  'F  =  fe  -  <y/  +  ^  .  ^  =  ry  /  c  is  the 

optical  wave  number,  ©  is  the  optical  wave's  angular  frequency  and  ^  is  the  optical  phase. 
Optical  fields  obey  Maxwell's  wave  equation: 


(23) 


where  is  the  total  transverse  beam  current.  If  the  expression  for  A ,  (22),  is  inserted 
into  (23)  and  the  left-hand  side  of  the  equation  is  written  out,  a  quite  long  and 
complicated  expression  results.  The  left-hand  side  can  be  simplified  greatly  by  making 
the  assumptions  that  the  optical  wave  varies  slowly  in  space  with  respect  to  an  optical 
wavelength  A ,  and  varies  slowly  in  time  with  respect  to  a  single  optical  period  =  A/ c . 
Mathematically,  this  means: 

dE  jr  d(f>  A  dE  T  r:  d(p 

—  «coE,—^«co<p  and  — kkKqE  , «kQ(f) .  (24) 

dt  dt  dz  dz 

Using  (24)  to  simplify,  (23)  reduces  to: 


2(dE 
c[  8t 


2E  f 


(25) 
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where  we  have  defined  the  two  orthogonal  unit  vectors  as: 

=  [cos  'F,-  sin  ^,0]  and  =  [siti  cos  'F,o] .  (26) 

Through  the  slowly-varying  amplitude  and  phase  approximation,  (25)  is  converted  from  a 
second-order  differential  equation  to  two  first-order  equations: 


2  ^ 


V^- 


1  d 
c'  at' 


dt 


—  and 

c 


^  2  1  2Ed(p  Ak  -  . 

^  = - 

c  dt  c 


V'  — 


c'  at' 


(27) 


(28) 


Now  attention  is  turned  to  the  right-hand  side  of  Maxwell's  equation.  The  current 
of  a  single  electron  is  J,  =  -ecfi^S^  (^-^)  where  e  is  the  electron  charge  magnitude, 
is  the  dimensionless  transverse  electron  velocity  (8),  J'(3c-i;.)  is  the  three 

dimensional  Dirac  fimction,  and  r,  is  the  position  of  the  i'th  electron  [4],  Expanded  out, 
the  transverse  current  for  one  electron  becomes: 


^  6^  B 

Jjj= - [cos(A:oZ),  sin(A:o  z),0](J^  (x-r^), 

ymco 


(29) 


where  B  is  the  magnitude  of  the  undulator's  magnetic  field,  y  is  the  Lorentz  factor  (3),  m 
is  the  electron  mass,  and  A^^is  the  undulator's  wave  number  (1).  Substituting  (29)  into 
(27)  and  (28),  the  wave  equation  becomes: 


—  D 

= - cos(T'  +  kQ z)S^ (x-rj) 

ymco 


J  •£  = 


e^B  . 


ymco 


sin('F  +  kQ  z)S^  (x—rj) 


(30) 


(31) 
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Equations  (30)  and  (31)  deal  only  with  a  single  electron.  The  total  transverse 
current  is  the  sum  over  all  electrons.  This  is  represented  by  averaging  (...)  over  a  fixed 
voliune  weighted  by  the  electron  density  p  resulting  in: 


ydt  ) 


ItjbK 

r 


p  {cos(^ +({)))  and 


(32) 


7jbK 

\c  dt )  y 


p(sinC'+^zJ)), 


(33) 


where  K  is  the  undulator  parameter  (10)  and  ^  is  the  electron  phase.  Equations  (32)  and 
(33)  can  be  simplified  by  defining  the  dimensionless  current: 


._mp{e7KL)^ 

J  ~  3  2  ’ 

y  me 


(34) 


where  L  is  the  undulator  length  and  N  is  the  number  of  undulator  periods  (21),  and  the 
dimensionless  optical  field  strength  |a|  in  (21).  The  results  are: 


a  =  -y  (cos(^  +  <f)ff  and 


<^  =  ^{sin(^  +  <^)), 


(35) 

(36) 


where  a  =  \a\e'^  •  Equations  (35)  and  (36)  can  be  combined  back  into  one  equation  in 
complex  form: 


(37) 
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C.  LOW  GAIN  DERIVATION 


Gain  is  the  fractional  change  in  power  of  the  optical  wave  as  it  makes  a  single 
pass  down  the  undulator.  As  mentioned  in  the  introduction,  there  are  two  basic  PEL 
designs:  oscillators  and  amplifiers.  Oscillators  typically  use  low  current  ( y  <  1 )  to  build 
the  optical  field  through  multiple  passes  of  the  light  wave  through  the  undulator  and 
between  the  mirrors  of  an  optical  resonator.  Amplifiers  on  the  other  hand  typically  rely 
upon  high  current  (  7  » 1 )  to  provide  high  gain  in  a  single  pass.  The  PEL  under 
construction  at  Thomas  Jefferson  National  Accelerator  Pacility  is  an  oscillator  design  so 
only  low  gain  is  discussed  here  [5].  See  reference  [6]  for  an  explanation  of  high  gain 
behavior. 

Using  conservation  of  energy,  any  change  in  the  optical  wave  energy  is  a  result  of 
an  opposite  change  in  the  electron  beam  energy  so  that  the  gain  can  be  calculated  by 
dividing  the  change  in  electron  beam  energy  by  the  optical  wave  energy.  To  find  this 
overall  change  in  electron  beam  energy  (which  consists  of  electrons:  some  gaining  and 
some  losing  energy  (21)),  we  must  first  find  an  expression  for  the  energy  change  of  a 

single  electron,  Aymc^ . 

To  help  understand  single  electron  behavior,  we  define  the  dimensionless  electron 
phase  velocity  y  as  the  dimensionless  time  r  derivative  of  the  electron  phase  4' ;  thus 

o 

(;  =  v  =  L[{ko  +  ku)fi2  ~^o\-  ^  electron's  change  in  energy  is  related  to  its  change  in 
phase  velocity  At> .  Since  k»ko,  Av  =  ■  Using  (13)  and  (17),  the  relationship 

between  Ay  and  Av  is: 
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f  X  \ 

Av  =  LkAjSz  =  NX^kAfiz  =  IttN  Afi^  *  2;rA^| 

\  X  ) 

The  energy  change  of  a  single  electron  is  then: 


^  2y^  ^ 
^1+a:  j 


Ap2«A7iN^.  (38) 
y 


Aymc  = 


2  ymc^Av 


AttN 


(39) 


In  the  case  where  the  optical  field  is  weak  {a  «n:)  and  the  gain  is  low  {j  «n: ), 
there  is  little  change  in  uover  the  entire  undulator  length.  Therefore,  using  perturbation 
theory,  a  solution  to  the  pendulum  equation  (21)  expanded  in  powers  of  the  initial  optical 
field  strength  : 


=  +'^o'^-^bos(Co  +i^o^')-cos^o+Vo'Z'sin^ol+"->  (40) 

and  the  electron  phase  velocity  becomes: 


u(t,C„  )  =  u„+  —  [sin(4'o  + 1^0  r)  -  sin  ]+ 


o. 


^2 


O  L 


-^{cos(2^iy  +2U(,r) -cos 2^(2 )+ cos Upr-l-u^r sin cos(^o  +^o‘^) 


(41) 


where  and  Oq  represent  initial  values  of  their  respective  variables. 

The  average  phase  velocity  of  all  sampled  electrons  within  the  beam  (pj  gives  the 

change  from  the  initial  phase  velocity,  Ao  =  (u)-Uq.  This  makes  the  average  change  in 
beam  energy: 


—  2 

Aymc  = - 


AttN 


(42) 


Equation  (41)  can  be  converted  by  averaging  over  the  initial  phases  : 
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=  o(T,Co)dCo=^-  COSVoT-l  +  l-UoTSmUoT  .  ( 

ITT  OqI  2 

We  must  also  include  the  fact  that  the  electron  beam  radius  does  not  equal  the 
optical  beam  waist.  Conservation  of  energy  only  applies  within  the  area  where  both 
beams  overlap.  Therefore,  we  must  introduce  a  filling  factor: 


where  is  the  electron  beam  radius  and  is  the  radius  of  the  optical  beam  waist.  For 
the  whole  beam,  the  number  of  electrons  is  the  electron  density  p  multiplied  by  the 


volume  dV occupied  by  the  optical  beam  one  wavelength  long  so  that  (42)  becomes: 


Aymc^  = 


pFdVymc^  ((u)  -Vq) 


The  optical  energy  within  the  volume  dV is  E^dV  I  Ait  so  the  low  gain  equation  is: 


A/mc^  _  f pFdVymc ^  ({o) -  Vq  ) V  4;r 


Energyop„cAL 


E^dV 


Inserting  (34)  and  (43)  into  (46),  the  gain  finally  becomes: 


j  \l-2cosvQT-VQrsinvQr\ 


where  Fhas  been  absorbed  into  j  for  brevity.  To  gain  a  better  imderstanding  of  equation 


(47),  Figure  5  is  a  plot  of  final  ( r  =  1 )  gain  versus  the  initial  electron  phase  Vq  with 7=1 . 
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G 


Figure  5.  Low  Gain  Spectrum  for  Low  Current  and  Weak  Fields 

Notice  that  the  gain  spectrum  is  anti-symmetric  about  resonance  =  0 .  Why 
does  resonance  result  in  zero  gain?  Electrons  enter  the  undulator  with  a  random  initial 
phase  from  -nil  to  'inH.  When  an  electron's  phase  is  7r/2  <  <  37i/2,  equation 

(21)  says  the  phase  aeceleration  is  negative.  Negative  phase  acceleration  makes  the 
change  in  phase  velocity  Al>  also  negative.  Equation  (39)  states  that  when  an  electron 
has  -  Au  ,  it  losses  energy.  That  energy  is  transferred  to  the  optical  field  thereby 
increasing  the  laser's  gain.  On  the  to  other  hand,  when  an  electron's  phase  is  -nil  < 

<  nil,  then  the  phase  acceleration  is  positive  so  the  electron  is  gaining  energy.  The  only 
way  for  the  electron  to  increase  in  energy  is  by  absorbing  energy  from  optical  wave 
thereby  decreasing  the  laser's  gain.  If  the  electrons'  initial  velocity  is  at  resonance,  they 
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bunch  together  dX  ^  12  and  equal  number  of  electrons  gain  and  lose  energy.  The  net 

effect  is  zero  optical  gain.  Figure  6  graphically  shows  the  resonance  condition. 

Plotted  in  the  main  window  is  the  phase  space  path  (^versus  l>)  of  30  electrons 
as  they  travel  down  the  undulator  (their  paths  darken  as  they  approach  r  =  1 ).  Those 
electrons  that  are  decreasing  in  ware  losing  energy  (increasing  the  laser's  optical  energy) 
and  those  that  are  increasing  in  v  are  gaining  energy  (depleting  the  laser's  optical 
energy).  The  numbers  of  electrons  gaining  and  losing  are  approximately  equal,  and  they 
are  all  approaching  the  same  phase  of  nil.  The  top  right  window  depicts  gain  growth 
as  a  function  of  r ,  and  it  never  rises  above  zero.  The  lower  right  graph  shows  the 
evolution  of  the  optical  phase  ^(r) . 

***  FBi;^  Phase  ■  Space  Evolution  *** 
j==l  ''o“^  N=40  ' 

0=0 

0.13 


0 

0.1 

0 

-x/2  ^  3x/2  Ox  1 


Figure  6.  Phase  Space  Diagram  of  the  Resonance  Condition 
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While  the  resonance  condition  is  a  good  case  in  understanding  electron  behavior 
and  bunching,  it  is  pointless  designing  a  laser  with  zero  gain.  On  the  other  hand,  (47)  has 
a  maximum  positive  gain  at  Vq  =  2.6 .  To  help  understand  why.  Figure  7  illustrates  this 
case.  Here,  more  electrons  lose  phase  velocity  than  gain,  and  instead  of  bunching  at  nil, 
they  bunch  at  n .  This  can  be  seen  analytically  by  combining  (20)  and  (39)  to  yield  the 
fractional  change  in  electron  energy: 


L 

7 


AnN 


COS^O' 


(48) 


Equation  (48)  has  the  most  negative  value  when  ^Q-n.  Consequently  tvith  j  =  l  in 
Figure  7,  ln(l  +  G)  peaks  to  a  maximum  value  of  0.13,  representing  13%  gain. 


**♦  FBL  Phase  Space  Evolution  *** 

j=l  a  =1  V  =2-5  N=40 

o  o 

0=0 


0.13 


0 


0.1 


0 


Figure  7.  Phase  Space  Diagram  of  Maximum  Gain  for  Low  Current  &  Weak  Fields 
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D.  GAIN  DEGRADATION  DUE  TO  ELECTRON  BEAM  QUALITY 

Until  now,  a  perfect  electron  beam  entering  the  undulator  has  been  assumed.  This 
means  that  all  electrons  enter  the  imdulator  with  the  same  initial  energy  (i.e.  the  energy 
spread  is  Ay  /  y  =  0 )  and  they  enter  with  zero  displacement  and  angular  error  from  the 
undulator's  z-axis  (i.e.  the  emittance  £■  =  0). 

Concerning  PEL  performance,  (47)  states  gain  increases  with  increasing  current 

density  j.  Using  (16)  and  (34),  j  oc  ,  where  /is  the  electron  current,  so  there  are 

three  ways  to  increase  j.  However,  A,  is  usually  fixed  for  a  particular  application.  In 
addition,  increasing  /  tends  to  decrease  the  beam  quality  from  the  accelerator,  and 
increasing  N  tends  to  increase  the  FEL's  sensitivity  to  beam  quality  by  narrowing  the  gain 
spectrum  bandwidth.  So,  there  exists  a  trade-off  between  high  current  density  and  good 
beam  quality  [7]. 

To  help  quantify  how  poor  beam  quality  affects  PEL  performance,  a  Gaussian 
distribution  is  assumed  for  the  energy  spread  using  (39)  to  define  the  standard  deviation 
as: 

<Tq  =  Ay  =  4;zi\/Ay  /  y ,  (49) 

If  two  electrons  enter  the  undulator  with  the  same  electron  phase  ^ ,  they  will  drift  apart 
by  A^  =  AuAt  .  If  cTp  =  Ay  « ;r ,  then  the  phase  separation  is  as  much  as  XI 2  by  the 
time  the  two  electrons  leave  the  undulator.  This  impairs  bunching  significantly  reducing 
gain.  The  typical  energy  spreads  from  a  good  accelerator  is  on  the  order  of  10"^ . 

Gaussians  are  also  assumed  to  represent  displacement  and  angular  errors.  Their 
associated  standard  deviations  are: 
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and  (50) 

where  is  the  electron's  rms  distance  from  the  z-axis  and  9  is  the  electron's  rms 

injection  angle.  As  with  cr^ ,  when  either  of  these  standard  deviations  is  near  n ,  gain 
begins  to  degrade. 


E.  SHORT  PULSE  EFFECTS 

Most  FEL  oscillators  use  short  electron  pulses  rather  than  a  continuous  beam. 
These  electron  pulses  in  turn  produce  a  short  optical  pulse.  At  the  undulator  entrance 
(7  =  0),  electron  and  optical  pulses  enter  together  where  the  electron  pulse  gives  up  a 
portion  of  its  energy.  Although  both  pulses  may  enter  the  undulator  at  the  exact  same 
time  (i.e.  synchronized),  their  relative  positions  change  as  they  travel  down  the  undulator. 
The  electron  pulse  travels  somewhat  slower  than  the  speed  of  light  so  it  slips  back 
relative  to  its  starting  point  with  the  optical  pulse.  This  slippage  distance  is  approximated 
by  assuming  the  FEL  is  at  resonance.  Resonance,  by  definition,  is  when  one  wavelength 
of  light  passes  over  an  electron  as  that  same  electron  travels  one  undulator  wavelength 
Xq  .  By  the  time  the  electron  exits  the  undulator,  it  has  slipped  back  a  distance  of  'NXq  . 

If  the  electron  pulse  length  is  comparable  to  this  distance,  then  the  electron-photon 
interaction  is  dominated  by  short  pulse  effects  [8]. 

If  the  distance  between  mirrors  is  5,  then  an  optical  pulse  arrives  at  the  undulator 
entrance  ( r  =  0 )  in  time  intervals  of  2Slc.  As  the  optical  pulse  enters  the  imdulator,  the 
distance  ahead  of  its  associated  electron  pulse  is  the  FEL's  desychronism  d.  It  is 


normalized  to  the  slippage  distance,  and  equals: 
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(51) 


2AS 

a  = - 

NXo 

Operationally,  d  is  changed  by  moving  one  of  the  mirrors  in^vard  by  a  small  amount  AS . 
At  first  thought,  d=0  (i.e.  synchronized  pulses)  would  seem  to  result  in  the  most  coupling 
since  the  peak  of  the  electron  pulse  would  coincide  with  the  optical  pulse  peak  each  pass. 
However,  this  is  not  the  case.  In  fact,  it  results  in  zero  final  power!  Figure  8  shows  the 
pulse  evolution  of  such  a  case  where  the  optical  power  is  evolving  to  zero. 


1'  se'-.-EvoiiiMon  *  * 

■Qrlp; 

•:4=Q 

,  ^4=^'PXO:O01. 

-1.5  z  1.5  -31  V  31  -31  V  31 


-1.5  z  1.5-31  V  310  n  300 


Figure  8.  Synchronized  Pulse  Evolution 


To  help  understand  how  the  electron  and  optical  pulses  evolve,  the  longitudinal 
coordinate  is  scaled  to  the  slippage  distance  by  z  /  NX  z .  This  means  the  optical 

pulse,  travelling  at  speed  c,  remains  fixed  in  z  while  the  electron  pulse,  travelling  slower 
at  ,  slips  back  in  z.  The  lower-left  picture  of  Figure  8  is  a  graph  of  the  dimensionless 
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current  density  j  versus  z  at  r  =  0  (dark  gray)  and  at  r  =  1  (light  gray).  The  electron 
pulse  is  assumed  to  be  parabolic  of  the  form  j{z)  =  {\-2z^  /  crl)]  for  |z|  <(7^1^  and 

zero  elsewhere.  In  Figure  8,  =  1  and  has  a  maximum  value  of  j=2. 

The  number  of  passes  through  the  undulator  is  represented  by  n=300  and  the 
cavity  power  loss  per  pass  is  VQ  =  1/20.  Noise  is  simulated  in  Figure  8  by  inserting  a 

random  phase  displacement  per  electron  with  rms  value  of  =  \0~* .  Just  above  the 
electron  pulse  graph  is  shown  the  dimensionless  optical  amplitude  a  evolution  as  a 
function  of  both  z  and  n.  The  bottom,  middle  picture  graphs  (47)  showing  gain  G  as  a 
function  of  electron  phase  velocity  o .  Just  above  it  is  the  FEL's  power  spectrum  as  a 
function  of  n.  The  middle  right  graph  it  is  the  o  distribution /as  a  function  of  n  while  the 
graph  just  below  it  depicts  dimensionless  power  P  as  a  function  of  n.  The  power  P  is  the 
square  of  the  dimensionless  optical  pulse  amplitude  a,  and  can  be  converted  into  a  real 
average  power  in  four  steps.  First,  using  the  definition  of  a  in  (21),  the  real  electric  field 
magnitude  E  is  found.  Second,  the  power  density  is  now  known  since  it  equals  l%7t . 
Third,  the  peak  power  is  the  power  density  multiplied  by  the  beam  area,  and  last,  the 
average  power  is  the  peak  power  multiplied  by  the  duty  cycle. 

The  optical  pulse  initially  grows  from  noise,  and  is  as  wide  as  plus  the 
slippage  distance  Az  =  1 .  Once  the  optical  pulse  form,  the  electron  pulse  couples  with  it, 
and  it  intensifies  (see  the  lower  right  P(n)  graph  in  Figure  8).  It  also  narrows  (see  the 

middle  left  plot  of  \a(z,  n)| ).  However,  since  the  electron  pulse  drifts  back,  it  amplifies 
the  optical  pulse's  trailing  edge  on  every  pass  through  the  undulator.  Equation  (47) 
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shows  this  analytically.  Very  little  gain  occurs  in  the  first  third  of  the  undulator  (r«l), 
and  the  majority  of  gain  occurs  towards  the  end  (  r  « 1 ).  As  stated  earlier,  gain  requires 
electron  bunching,  and  bunching  takes  time.  The  graph  of  \a{z,  «)|  shows  this  as  the 
optical  pulse  peak  decreases  in  z  with  each  increase  in  n.  Eventually,  the  optical  pulse 
moves  back  so  far  that  it  decouples  from  the  incoming  electron  pulses,  and  dies  out 
aroimd  n=300. 

To  ensure  the  optical  and  electron  pulses  remain  coupled,  the  optical  pulse 
pathlength  As  must  be  shortened  (i.e.  J  >  0 ).  Even  a  small  desynchronism  can 
accomplish  the  desire  effect.  Figure  9  uses  the  same  parameters  as  those  in  Figure  8  with 
the  desynchronism  of  d=0.003.  Because  the  optical  pulse  is  advanced  for  each  pass,  its 
position  in  z  remains  constant.  This  allows  the  electron  pulse  to  "sweep"  across  the  entire 
optical  pulse  leaving  it  undistorted.  As  a  result,  the  FEL's  power  grows  and  eventually 
levels  out  to  a  steady-state  value. 

Unfortunately,  increasing  ci  indefinitely  is  coimter-productive.  With  increasing  d, 
eventually  the  optical  pulse  is  advanced  so  far  that  it  again  decouples  with  the  incoming 
electron  pulse  reducing  steady-state  P.  For  large  desynchronism,  the  optical  field  has  a 
long  exponential  tail  of  the  form  [8]: 

|a(z)|oce-"''2^  (52) 
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FBL  Pulse  EJvoiTxtion  **** 

3=2 

xr^=l 

0^=0  0^=6  d=0-003 

Q=20 

A=0 

A=0  8=0  N=24  8^=0.0001 

-1.5  z  1-5  -31  V  31  0  n  300 


Figure  9.  Desynchronized  Pulse  Evolution 


F.  THE  KLYSTRON  UNDULATOR 


The  klystron  is  a  two-stage  undulator  specifically  designed  to  improve  FEL  gain 
in  weak  optical  fields  [9].  Figure  10  shows  the  design  of  a  klystron  undulator. 

Simply  stated,  it  is  nothing  more  than  an  conventional  undulator  cut  in  half  and  spread 
apart  to  create  a  drift  space  between  sections.  The  first  section,  called  the  "modulator". 


Entrance 


Exit 


Modulator 


Drift  Space  Radiator 


J 


Figure  10.  Kylstron  Undulator 
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prepares  the  electrons  for  bunching  by  developing  phase  velocity  differences.  The 
middle  section  between  the  conductors  provides  space  for  the  electrons  to  bunch  as  they 
drift.  The  final  section  is  called  the  "radiator"  where  highly  bunched  electrons  create  and 
amplify  the  optical  field.  The  drift  section  can  be  shortened  by  using  bending  magnets  (a 
dispersive  section)  to  reduce  overall  length  [9]. 

The  drift  or  dispersive  section  strength  is  given  by  the  dimensionless  parameter: 

N 

(53) 

where  N  is  the  number  of  undulator  periods  in  both  the  modulator  and  radiator,  and  is 
the  drift  space  length  in  number  of  undulator  periods.  Within  the  drift  space,  o  is 
constant  so  the  only  electron  phase  change  is  =  uD .  Electron  bunching  will  therefore 
occur  without  optical  pulse  interaction  so  the  optical  field  remains  constant. 

Via  simulation,  the  low  current  gain  curve  as  a  function  of  Vq  can  be  constructed. 
Figure  1 1  shows  the  gain  spectrum  with  j=l  and  D=1 .  Comparing  with  Figxire  5,  the  gain 
increases  from  0.13  to  0.37,  about  three  times  as  much.  For  jLxl,  the  gain  equation  can 
be  approximated  near  resonance  as: 

G{Vo)^^siniVoD)  (54) 

The  maximum  gain  therefore  becomes  jDIA  when  Ug  =  tc!  2D . 
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Gain  and  Phase  Curves  **** 
j=l  a^=l  <r=0 

D=1  N=100  6^=0 


Figure  11.  Low  Current  Klystron  Gain  Versus  Initial  Phase  Velocity 


A  negative  aspect  to  using  a  klystron  is  that  the  electron  beam  quality  must  be 
sufficient  to  keep  AVq  <7tl  D  .  Otherwise,  the  gain  is  negative.  This  is  a  tighter 
constraint  than  that  placed  upon  conventional  undulators  (  <7t),  and  can  become  the 

limiting  design  factor  when  D  is  large. 
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III.  TJNAF  FREE  ELECTRON  LASER 


A.  KWFEL  DESIGN 

As  mentioned  in  the  introduction,  TJNAF  is  planning  to  achieve  first  light  from 
their  kW  FEL  this  summer.  Figure  12  is  a  basic  diagram  of  this  laser.  Free  electrons  are 


created  from  the  350  KeV  photocathode  gun  and  initially  accelerated  in  the  10  MeV 
cryomodule.  At  this  point,  they  are  accelerated  further  with  the  32  MeV  cryomodule.  As 
shown  in  the  theory  section,  shortening  the  pulse  width  increases  j  that  in  turn  improves 
the  gain  G.  To  capitalize  on  this,  electrons  are  "pre-bunched"  using  a  compressor  before 
entering  the  undulator.  Upon  leaving  the  imdulator,  electrons  are  "unbunched"  using  a 
decompressor  before  they  enter  the  first  major  bend.  This  helps  reduce  energy  loss 
whenever  electrons  change  direction.  Electrons  are  eventually  sent  to  one  of  the  three 
dumps  dependant  upon  the  system  configuration.  The  two  10  MeV  dumps  are  used  when 
testing  the  electron  gun  and  cyromodules.  During  normal  operation,  the  electrons  are 
sent  to  the  42  MeV  dump  following  the  decompressor.  The  recirculation  ring  serves  an 
important  purpose  by  redirecting  the  electrons  back  through  the  32  MeV  cryomodule  out 
of  phase.  This  serves  two  purposes.  First,  it  conserves  energy  by  using  the  recirculated 
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electrons  to  power  in  part  the  cryomodule.  This  reduces  what  the  cryomodule  would 
otherwise  need  in  RF  power.  Second,  it  reduces  the  beam  dump  size  by  removing  32 
MeV  from  the  electrons  before  discarding  them. 


B.  PROPOSED  20  KW  PEL  PARAMETERS 

Following  successful  demonstration  of  their  kW  FEL,  TJNAF  plans  to  improve 
their  design  to  provide  20  kW  on  target.  Table  1  below  is  a  list  of  proposed  parameters  to 
meet  this  goal,  and  was  used  in  the  simulation  work  contained  herein. 


PARAMETER 

Optical  Wavelength  A 

1.05  pm 

Undulator  Wavelength 

20  cm 

Peak  Undulator  Magnetic  Field  B 

2kG 

Undulator  Periods  N 

24 

Undulator  Length  L 

6  m 

Undulator  Parameter  K 

0.79 

Dispersive  Strength  D 

1,2  or  3 

Electron  Energy  ymc^ 

200  MeV 

Maximum  allowed  Energy  Spread  Ay  /  y 

6% 

Electron  Beam  Radius 

0.5  mm 

Electron  Pulse  Length  cr^ 

2.0  ps 

Initial  Phase  Velocity  Vq 

0.0  (on  resonance) 

Initial  Electron  Phase  Velocity 

Standard  Deviation  cr^ 

0.3 

Electron  Injection  Angle  Standard  Deviation  Cg 

0.15 

Peak  Current  I 

200  A 

Average  Current 

8  mA 

Current  Density  j 

2 

Initial  Optical  Field 

0.0 

Cavity  losses  MQ 

0.1 

Peak  Optical  Power  within  cavity 

40  GW 

Repetition  Rate 

74.85  MHz 

Duty  Cycle 

4jdfr^ 

Average  Optical  Power  P  within  cavity 

1.6  MW 

Required  Efficiency  for  20  kW  on  target 

1  % 

Table  1.  TJNAF  20  kW  Proposed  Parameters 
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IV.  SIMULATION  RESULTS 


A.  KLYSTRON  GAIN 

The  first  design  goal  was  to  discover  the  steady-state  gain  behavior  as  a  function 
of  both  dispersive  strength  D  and  desynchronism  d  for  a  proposed  20  kW  klystron. 
TJNAF  provided  necessary  design  parameters  as  provided  in  Table  1 .  TJNAF  is 
considering  D's  of  1, 2  or  3.  LT  Richard  Steele  conducted  simulation  work  for  the  cases 
of  D=Q  (a  conventional  undulator,  no  klystron)  and  1  [10]  while  this  thesis  deals  with  the 
cases  D=2  and  3.  For  completeness,  the  final  results  of  both  are  combined  in  the  follow- 
on  graph.  Desynchronism  was  examined  fi:om  d=0.0  (i.e.  synchronized  pulses)  to  d=0.3. 

The  steady  state  gain  behavior  was  revealed  using  the  pulse  evolution  output  as 
described  in  Chapter  II,  Section  E  with  one  exception.  Instead  of  plotting  power  P  as  a 
function  of  the  number  of  passes  n  down  the  undulator,  P  is  replaced  with  gain  G  as  a 
function  of  n.  The  simulation  was  run  only  long  enough  for  gain  to  reach  a  steady  state 
value.  The  results  are  shown  in  Figure  13.  Desynchronism  is  plotted  along  the  x-axis 
with  G  as  the  y-axis.  Simulations  were  conducted  in  steps  of  d=0.02  for  all  four 
dispersive  strength  cases.  The  only  deviation  fi’om  Table  1  in  design  parameters 
concerned  the  mirror  losses  \/Q.  To  shorten  the  time  required  for  each  simulation,  mirror 
losses  where  neglected  by  making  ^lxl0'°.  This  allowed  G  to  be  evaluated  accurately. 
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steady  State  Gain 


Figure  13.  Steady  State  Gain  For  Various  Dispersive  Strengths 

For  the  conventional  undulator  (D=0)  starting  at  synchronized  pulses  (d=0),  the 
steady-state  gain  begins  at  zero,  and  then  increases  to  a  maximum  value  of  G=0.18  at 
d=0.03.  From  there,  any  further  increases  in  reduce  peak  gain  imtil  it  again  is  zero  at 
d=0.l3. 

As  shown  in  Chapter  II,  Section  F,  a  klystron  improves  gain  in  the  presence  of 
weak  optical  fields.  Figure  13  quantifies  this  improvement.  With  D=l,  the  klystron 
provides  a  2.5  increase  in  peak  gain  (from  G=0. 1 8  to  G=0.50).  One  unfortunate 
consequence  is  that  maximum  gain  now  occurs  at  a  much  larger  value  of  d=0. 1 1 . 
Increasing  c?  decreases  final  power  further,  as  will  be  seen  later.  A  dispersive  strength  of 
D=2  actually  provides  the  most  gain;  G=0.67.  Notice  how  the  maximum  gain  now 
occurs  at  the  same  d  for  each  D  value  of  1, 2  or  3. 
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An  interesting  phenomenon  occurs  from  D=2  to  D=3>.  Although  dispersive 
strength  increases,  there  is  no  further  increase  in  steady  state  gain.  Therefore,  there 
appears  to  be  an  estimated  optimum  value  of  Z)=2. 

B.  KLYSTRON  POWER  AND  ENERGY  SPREAD 

The  second  design  goal  is  to  define  the  final  average  optical  power  as  functions  of 
both  d  and  D.  Table  1  provided  the  operating  parameters,  and  the  same  values  as  used  in 
gain  for  Z)=0,  1, 2  and  3,  d  in  step  size  of  0.02  between  <^=0.0  to  <#=0.3.  Again,  LT 
Richard  Steele  ran  simulations  for  values  D=0  and  1  [10],  and  his  final  results  are 
included  in  the  following  results. 

The  simulated  optical  pulse  was  allowed  to  evolve  to  steady  state  requiring  n= 
2,000  passes  for  small  d  to  «=300  passes  for  large  d.  The  final  dimensionless  power  was 
read  from  each  output  graph,  and  then  converted  into  a  real  average  power.  The  residts 
are  shown  in  Figure  14.  Desynchronism  remains  on  the  x-axis  with  average  power  on  the 
y-axis. 

From  Figure  14,  a  conventional  undulator  with  Z)=0  provides  the  most  power. 
Maximum  power  occurs  at  a  small  value  of  <#=0.01  (see  Chapter  II,  Section  E),  and  the 
conventional  undulator  provides  25  times  as  much  as  the  klystron  case  for  D=\ .  As 
mentioned  earlier,  a  klystron  improves  gain  in  weak  fields,  but  imfortunately,  it  is  at  the 
expense  of  final  power  at  saturation. 
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Final  Average  Power 


Figure  14.  Final  Average  Power  For  Various  Dispersive  Strengths 


Klystron  Final  Average  Power 


Figure  15.  Klystron  Final  Average  Power  For  Various  Dispersive  Strengths 
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Since  the  conventional  undulator  provides  much  more  power,  why  use  a  klystron 
at  all?  Higher  powered  lasers  (kW  and  above)  may  use  higher  cavity  to  increase  output. 
Therefore,  the  gain  G  must  be  improved  to  overcome  these  losses.  The  first  way  to 
increase  G  is  to  increase  j,  but  that  is  limited  by  the  electron  gun  technology.  The 
klystron  provides  a  secondary  means  to  provide  the  additional  gain  required.  So  while 
the  klystron  lowers  output  power,  it  may  very  well  be  a  necessary  evil  in  achieving  the 
required  level  of  gain  to  make  the  laser  work  in  the  first  place. 

Since  the  conventional  imdulator  so  dominates  the  graph.  Figure  15  eliminates 
this  case,  and  allows  comparison  in  terms  of  power  of  differing  dispersive  strengths  D. 
Here,  D=l  provides  about  twice  as  much  power  as  the  case  D=2  with  diminishing  returns 
shown  for  J>=3.  While  the  cases  of  D=0, 2  and  3  show  an  early  maximum  then  tailing  off 
to  zero,  £>=1  exhibits  quite  irregular  behavior.  After  peaking  at  <£=0.01,  it  decreases  until 
d=0. 1 3  where  then  it  improves  to  its  highest  value  at  d=0. 1 7.  The  reason  for  this  is 
explained  at  the  end  of  this  section. 

To  understand  a  little  about  the  power  behavior,  let  us  examine  a  few  points  on 
the  graph  in  detail  for  D=2.  First,  let  us  look  at  d=0.01  where  power  is  maxirmim. 

Figure  16  shows  the  pulse  evolution  for  this  point  on  the  graph. 
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FEIr  Pulse  Evolution  ♦*** 

3=2 

00=0. 15 0,2=0. 3  d=0.01 

Q=10 

A=0 

A=0  D=2  N=24  8^=0 

Figure  16.  Pulse  Evolution  For  if=0.01 


For  a  general  explanation  of  each  plot  within  Figure  16,  refer  to  Chapter  II, 
Section  E.  Notice  how  the  optical  pulse  evolves.  It  begins  with  one  peak,  but  when  n  is 
approximately  650,  a  second  peak  forms  with  equal  intensity  to  the  first  peak.  The  power 
reaches  steady-state  and  the  electron  v  spread  is  At>  =  6.5 . 

The  third  design  goal  is  maintaining  energy  spread  Ay  /  y  below  6%.  An  energy 
spread  of  6%  or  less  is  necessary  for  TJANF  to  successfully  recirculate  electrons  (see 
Chapter  III  for  the  benefits  of  recirculation).  Using  (39),  this  translates  into  an  energy 
spread  of  A/  /  y  =  2%  that  is  well  within  the  design  requirement  of  6%.  Surprisingly  for 
Z)=2  and  £>=3,  the  energy  spread  is  constant  at  2%  for  all  values  of  desynchronism.  On 
the  other  end  of  the  desynchronism  scale.  Figure  17  shows  the  pulse  evolution  for 
ff=0.13. 
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Pulse  Evolution  ♦*** 
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Figure  17.  Pulse  Evolution  For  rf=0.13 

Here,  there  is  only  one  peak  with  a  long  exponential  ramp  formed  ahead  of  the 
main  pulse.  While  the  power  P  reaches  steady-state  quite  early  at  «  « 100  passes,  it  is 
only  a  third  of  its  maximum  value  in  Figure  16  at  £^=0.01.  The  optical  pulse  retains  this 
kind  of  shape  with  further  increases  in  d,  but  declines  gradually  to  zero  amplitude  at 
ch03. 

In  between  the  d  values  of  0.01  and  0.13,  the  optical  pulse  exhibits  erratic 
performance.  Figure  18  shows  such  an  intermediate  case  with  d=Q.01.  The 
transformation  is  taking  place  from  a  two  peak  optical  pulse  to  a  one  peak  pulse.  This  is 
also  an  example  of  trapped  particle  instability  [1 1]. 
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Figure  18.  Pulse  Evolution  for  </=0.07 


The  trapped-particle  instability  occurs  when  electrons,  in  the  presence  of  strong 
optical  fields  ( a  >  ;r ),  become  trapped  within  potential  wells  in  phase  space.  These 

electrons  begin  to  oscillate  at  the  "synchrotron"  frequency  [11].  This  behavior 

causes  the  beam  current  to  oscillate  that  in  turn  drives  modes  at  sideband  frequencies  of 
Vq±Vs,  shifted  away  from  the  fundamental  wavelength  by  AX/  A  =  0^  / XnN .  These 

sidebands  can  be  seen  in  Figure  18’s  power  spectrum  P{v,n) . 

Now  recall  the  power  curve  for  the  case  D=\  in  Figure  15.  It  reaches  a  peak  at 
cNO.Ol,  declines,  and  then  increases  again  from  chOA  to  0.2  before  finally  diminishing  to 
zero  at  d=0.n.  From  [10],  examining  the  optical  pulse  shape  with  increasing  d  yields  an 
explanation  for  this  behavior.  Although  the  optical  pulse  amplitude  peak  decreases  in  the 
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region  <^=0. 1  to  0.2,  the  exponential  ramp  of  the  form  shown  in  (52)  is  extending  at  a 
greater  rate  as  d  increases.  This  leads  to  an  overall  increase  in  final  average  power. 
However,  increasing  c?  beyond  0.2  provides  a  rapid  decrease  in  peak  amplitude  too  great 
for  the  extended  pulse  shape  to  overcome.  A  rapid  decrease  in  average  power  results 
until  at  large  d  there  is  no  laser  power. 
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V.  CONCLUSIONS 


The  purpose  of  this  study  was  to  quantify  the  behaviors  of  steady-state  gain,  final 
average  power  and  energy  spread  as  functions  of  desynchronism  d  and  dispersive 
strength  D  for  a  proposed  20  kW  klystron  free  electron  laser  (PEL).  This  proposal  is 
under  consideration  by  the  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  to 
improve  the  kW  FEL  they  are  now  constructing.  Design  parameters  are  listed  in  Table  1 . 
These  conclusions  combine  the  work  contained  herein  with  reference  [10]. 

The  two  major  design  goals  were  an  average  optical  power  of  20  kW  while 
maintaining  a  final  energy  spread  of  6%  or  less.  The  conventional  undulator  with  a 
desynchronism  ch0.03  will  achieve  these  goals. 

Gain  peaks  in  desynchronism  at  <#=0.05  for  £>=0  and  at  <#=0.1 1  for  D=l,  2  and  3. 
Gain  is  proportional  to  dispersive  strength  up  to  an  estimated  optimal  value  of  D=2. 
Further  increases  in  D  result  in  diminishing  gain.  A  klystron  with  D=2  provides  30% 
more  gain  than  Z>=1  and  70%  more  than  a  conventional  undulator. 

The  conventional  undulator  provides  25  times  as  much  power  as  a  klystron  with 
£>=1 .  Overall,  power  is  inversely  proportional  to  D  with  D=\  providing  twice  as  much  as 
D=2  and  almost  4  times  as  much  as  D=3 .  For  D=\ ,  peak  power  occurs  at  <#=0. 1 7,  and 
falls  off  rapidly  to  zero  at  fif=0.25.  For  D=2  or  3,  the  undulator's  power  peaks  early  at 
J=0.01 ,  and  then  gradually  declines  to  zero  at  <#=0.3. 

Energy  spread  is  2%  for  D=2  or  3  for  all  values  of  desynchronism.  The  energy 
spread  exceeds  the  design  goal  of  6%  only  for  the  case  Z)=0  within  the  range  ^#=0.01  to 
c/=0.03. 
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